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Abstract. The equation of a motion of curves in the projective plane is deduced. Local flows 
are defined in terms of polynomial differential functions. A family of local flows inducing 
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is discussed. Local motions defined by the traveling wave cnoidal solutions of the fifth-order 
Kaup-Kupershmidt equation are described. 
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1 Introduction 

The interrelations between hierarchies of integrable non linear evolution equations and motions 
of curves have been widely investigated in the last decades, both in geometry and mathematical 
physics. In the seminal papers [12, 13, 29], Goldstein, Petrich and Nakayama, Segur, Wadati, 
showed that the mKdV hierarchy can be deduced from local motions of curves in the Euclidean 
plane. Later, this result was extended to other 2-dimensionaI geometries [4, 5, 6, 7, 8, 32, 33] or 
to higher-dimensional homogeneous spaces [1, 2, 15, 16, 19, 21, 27]. The invariant curve flows 
related to integrable hierarchies are induced by infinite-dimensional Hamiltonian systems defined 
by invariant functionals and geometric Poisson brackets on the space of differential invariants of 
parameterized curves [22, 23]. Another feature is the existence of finite-dimensional reductions 
leading to Liouville-integrable Hamiltonian systems. Typically, these reductions correspond to 
curves which evolve by congruences of the ambient space. They have both a variational and 
a Hamiltonian description: as extremals of a geometric variational problem defined by the 
conserved densities of the hierarchy, and as solutions of a finite-dimensional integrable contact 
Hamiltonian systems. Examples include local motions of curves in two-dimensional Riemannian 
space-forms [25], local motions of star-shaped curves in centro-affine geometry [26, 32] and local 
motions of null curves in 3-dimensional pseudo- Riemannian space forms [27]. In [4, 5, 6, 7], 
K.-S. Chu and C. Qu gave a rather complete account of the integrable hierarchies originated by 
local motions of curves in 2-dimensional Klein geometries. In particular, they showed that the 
fifth-order Kaup-Kupershmidt equation is induced by a local motion in centro-affine geometry 
and that modified versions of the fifth- and seventh-order Kaup-Kupershmidt equations can 
be related to local motions of curves in the projective plane. Our goal is to demonstrate the 
existence of a sequence of local motions of curves in projective plane inducing the entire Kaup- 
Kupershmidt hierarchy. We analyze the congruence curves of the flows and we investigate in 
more details the congruence curves associated to the cnoidal traveling wave solutions [41] of the 

*This paper is a contribution to the Special Issue "Symmetries of Differential Equations: Frames, Invariants 
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fifth-order Kaup-Kupershmidt equation. In particular, we show that the critical curves of the 
projective invariant functional [3, 28] are congruence curves of the first flow of the hierarchy. 

The material is organized into three sections. In the first section we collect basic facts 
about the Kaup-Kupershmidt hierarchy from the existing literature [10, 11, 18, 24, 34, 37, 39] 
and we discuss an alternative form of the equations of the hierarchy. We examine the cnoidal 
traveling wave solutions of the fifth-order Kaup-Kupershmidt equation [41] and we exhibit new 
traveling wave solutions in terms of Weierstrass elliptic functions. In the second section we 
recall the construction of the projective frame along a plane curve without inflection or sextatic 
points and we introduce the projective line element and the projective curvature [3, 14, 30, 40]. 
Subsequently, we use the Cartan's canonical frame [3, 28] to study the equations of a motion of 
curves in the projective plane. Consequently, we construct local motions in terms of differential 
polynomial functions and we deduce the existence of a sequence of local flows inducing the 
Kaup-Kupershmidt hierarchy. In the third section we focus on congruence motions of the flows. 
In particular, we explicitly implement the general process of integration to determine the motions 
of the critical curves of the projective invariant functional [3, 28]. The symbolical and numerical 
computations as well the graphics have been worked out with the software Mathematica 8. We 
adopt [20] as a reference for elliptic functions and integrals. 



2 The Kaup-Kupershmidt hierarchy 
2.1 Preliminaries and notations 

Let J(M, M) be the jet space of smooth functions n : M — )• M, equipped with the usual coordi- 
nates (s, M(o)) ^(1)) • • • ) '^{h)-, • • • )• The prolongation of a smooth function u is denoted by j{u). 
Similarly, if u{s,t) is a function of the variables s and t, its partial prolongation with respect 
to the s- variable will be denoted by js{u). A map p : J(M, M) — t- M is said to be a polynomial 
differential function if 

p(u) = P(n(o), «(!),...,«(,,)) Vug J(M,M), 

where P is a polynomial va. h + 1 variables. The algebra of polynomial differential functions, 
J[u], is endowed with the total derivative 

oo r^ 

X - op 

and the Euler operator 

For each p G Ker(iS) there is a unique D^^p G J[u] such that 

p = Z)(Z)-i(p)), Z)^^p|o = 0. 
We consider the linear subspace 

P[u] = {p G J[u] : £:(u(o)D=^p + Au%)Dp) = O}. (1) 

The image of P[u] by the total derivative is denoted by P'[u] C J[u]. Next we consider the 
integro-differential operators 



A(0, u) = + 2u(l)s + UsCp, 
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E{(p, u) = (j)3s + 8u(ps + 7n^0 + 2{u2s + 4.u^) / Hr + 2 {u(p2s + 'iu^(p)dr, 

Jo Jo 

Us P 1 

y Jo y 

/3 16 o 71 2 41 13 \ , /59 35 \ , 

/ 2 49 \ 2 1 

and the linear operators 

V: J[u]^J[u], J: P'[u]^J[u], 5: P[u] ^ J[u] 
defined by 

P(tD) = L'^ro + 2u^q)Dk) + u(i)rD, 

J(q) = D^q + 8u^o)Dq + 7u(i)q + 2(n(i) + 4uJo))D-iq + 20'^ {u^o)D\ + 4n2o)q), 

'^(p) = ^^/(i)D~^(u(o)-D^p + 4uJo)L'p) + ^(20uJo)n(i) + 25n(i)n(2) + 10u(o)n(3) +U(5))p 

/3 16 o 71 2 41 13 \ ^ /59 35 \ ^2 

+ V2 ^ ^""("^ ^ 18"(^) ^ y ^(0)^(2) + jDp+ l^y n(o)n(i) + — n(3) j D p 

+ (2nJo) + Y^^(2)) ^^P - 2n(i)L>^p + ^n(o)I)^p + ^I?^p. 
From the definition of the operators is clear that 

1^i^)\j{u) = A(tt)|j-(„),M), 5(p)|j(„) = e(p|j(„),n), J'(q)|j(„) = H(q|j(„),u), 
for every tt) G j[u], p G P[u], q G P'[u] and every u G C°°(M,M). 
Lemma 1. The operators D, J and S satisfy 

VJ{p) = 18S{D-^p) - 27p Vp G P'[u]. 
Proof. A direct computation shows that 

A{E{(f), u),u) = 180 (^J^ (pdr,u^ -27(1) Vn, (/> G C°°(M, M). 

This implies 

1^J{p)\j{u) = A{E{p\j^^),u),u) = ISep-i (p) !,-(„), ^x) - 27p|,-(,) 
= 185(Z)-i(p))|,(„)-27p|,(,), 
for every p G P[u] and every u G C°^(M,M). This yields the required result. ■ 

2.2 Construction of the hierarchy 

According to [11, 34, 39] there are three sequences 

{f)n}n6N Clm(£:), {qnjneN C P[u] , {pnjneN C J[u] 

of polynomial differential functions defined by the recursion formulae 

i)n+2 = JV{i)n), VH)n)=D{qn), [)n = ^(pn) (2) 

and by the initial data 

f)0 = 1, f)i = M(2) + 4U(o). 
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Definition 1. The Kaup-Kupershmidt hierarchy {ICn} is the sequence of evolution equations 
defined by 

ICn : Ut +'D[)n\j{u) = 0. 

In view of (2), the equations of the hierarchy can be written either in the Hamiltonian form 

ICn : Ut+'D£{pn)\j(u) = 0, 
or else in the conservation form 
/C„ : ut + Dqn = 0. 

Remark 1. The polynomial differential functions f)„, q„, pn and the equations of the hierarchy 
can be computed with any software of symbolic calculus (see Appendix A.l). For n = 1,2 we 
find 



f)i(u 

i)2(u 

qi(u 

q2(u 
Pi(u 

P2(U 



U{2) + 4u- 



(0)' 



32 



12n(2)U(o) + 6u(^) + U(4) + — M(o), 



20 



15 



«(0) + ^^(1) + 10^(0)^(2) + ^i(4), 



56 



49 



""(0) + 70n(o)n(^) + 56n(Q)U(2) + -^U(2) + 35'U(i)n(3) + 14u(o)n(4) 



i(6), 



^-^(0)2^(2) + -'U(o), 



^^(o)''^(4) + 4u(o)n(2) + 2u(o)n(i) + -uf^y 



Consequently, the first two equations of the hierarchy are 

Ut + lOuUSs + 25UsU2s + 20u^Us + U5s = 0, 



Ut + U7s + 14uU5s + 49UsU4:S + 84:U2sU3s + 252uUsU2s + 70Ug + 56u U3s + 

Definition 2. Denoting by [r] the integer part of r, we set 



224 



-n'^u, = 0. 



1 



(l + (-l)"), A, 



-(l + (-l)"))(-27)[t] 



and we define {t)„}„gN C P[u] by 



0„ = 18|^(-27)W 



2/i, 



n = 1, . . . , 



/i=0 



where XOq = 0, tDi = 1/2 and rD„ = qn-2, n > 1. 

Proposition 1. The equations of the hierarchy can be written in the form 
ICn : dtU + cS(t)„)|j^(„) + XnUs = 0. 



(3) 



Proof. For n = 1,2 the proposition can be checked by a direct computation. We prove (3) 
when n is odd. By induction, suppose that (3) is true for n = 2p — 1. Note that 



A 



2p-l 



A 



2p+l 



0, 



2p+l — i^2p-l 



+ 1. 
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By the inductive hypothesis we have 

K,2p^l{u) = Ut + D{(\2p-l)\j^(u) = Ut+S{'02p~l)\j,{u) =Ut + e(02p-l|j4„),^x), 

which implies 

D{(\2p-l)\j,{u) = 0(D2p-l|j4«),«)- 

Using Lemma 1 we find 

K,2p+i{u) = ut + VJ{Dc\2p-i)\j,{u) = ut + 185(q2p-i)|j4„) ~ 27Dq2p-i|j4n) 
= ut + 18G(q2p-i|j4«),u) - 27G(t)2p-i|j4„),n) 
= ut + e((18tt)2p+i - 27t)2p-i)|j4„), n). 

Using 



/ hp-i 
18n32p+i - 27D2p_i = 18 tt)2p+i - 27 ^ (-27)'^ro2p_i-2h, 

\ h=0 

fhp+i-1 
W2P+1- (-27)''+'tt)2p-l-2(h+l) 



h=0 

2p+l 



h=0 

we obtain 



= 18 tt)2p+l - Y (-27)''tt)2p+l-2h 
\ h=l 

l2p+l 

= 18 ^ ( — 27)^rD2p+l-2/i = f2p+l 



IC2p+l{u) = Ut + Q{\}2p+l\j,{u),u) = Ut + S{\32p+l)\j,{u) + Mp+lUg- 

Next we prove (3) when n is even. By induction, suppose that (3) is true for n = 2p. Note that 

A2P+2 = i-27)P+' = -27A2p, ^2p+2 =p = i2p + l. 

By the inductive hypothesis we have 

IC2p{u) =Ut + D(q2p)|j4„) = Ut +5(t)2p)|j4u) - A2pMs = Ut + 0(t)2p|j4„),'w) + A2p'Us, 

which implies 

Diq2p)\js{u) = 'S>{^2p\j,{u},u) + \2pUs- 

From Lemma 1 we have 

^2p+2(^i) = ut + VJ{Dc\2p)\j,{u) =ut + 185(q2p)|j4„) ~ 27Z:>q2p|j 
= ut + 18G(q2p|j4„),'u) - 27(e(02p|j4„),^^) + ^sp^^s) 
= ut + 9((18tt)2p+2 - 27t)2p)|j4„), n) + -^2p+2Ms- 

Using 



18tD2p+2 - 27t)2p = 18 tD2p+2 - 27 ^(-27)''ra2p-2/i 
V h=0 



6 



E. Musso 



= 18 tt)2p+2 - (-27)''+^tt)2p-2h =18 tt)2p+2 - {-'^^f^2p+2-2h 

\ h=0 J \ h=l 

= 18 Yj (~27)'^rD2p+2-2fe = t)2p+2 

we find 

K.2p+2{u) = Ut + Q{V2p+l\j,{u),u)) + X2p+2Us = Ut + 5(t)2p+2) | {«) + >>2p+2Us- ■ 

2.3 Traveling waves of the fifth-order Kaup— Kupershmidt equation 

Several classes of traveling wave solutions of the fifth-order equation fCi have been considered in 
the literature [17, 38, 41]. In this section we generalize the elliptic families examined in [41]. The 
hyperbolic traveling waves found in [38] can be obtained as limiting cases, when the parameter 
of the elliptic functions tends to 1. First consider the third-order ODE 

k'" + 8kk' = 0, (4) 

where k', k" and k'" denote the first-, second- and third-order derivatives of a real-valued func- 
tion k with respect to the independent variable. The same notation will be used for vector- valued 
functions. Integrating twice (4) we find 

k'' = -\k' + l92k-\g,, (5) 
where g2 and 53 are real constants. Every k satisfying (5) generates the traveling wave solution 

u{s,t) = k(s - ^g2t 



of the first equation of the hierarchy. Clearly, (5) can be integrated in terms of the Weierstrass p 
functions, namely: if A{g2,g3) = — + 27(7| > 0, then 

3 

k{s) = — -p(s + c), s £ (2na;i — c, (2n + l)uJi — c), n £ Z, 

where cji is the real half period and c is a real constant. If A{g2,g3) < 0, then there are two 
types of solutions: 

3 

k{s) = —-p{s + c), s £ (2na;i — c, (2n + l)uji — c), n £ Z, 
3 

= --p(s + W3 + c), s£R, 

where wi and ^3 are the real and the purely imaginary half periods. When A(g2,<?3) < 0, the 
Weierstrass functions can be written in terms of Jacobi elliptic functions and we get 



13 1 

k{s) = -(1 + m) — -ns(s + c|m)^, k{s) = - (l — 2m + 3mcn(s + c|m)^) . 



The parameter m £ (0, 1) is 63 — 62, where ei > 62 > 63 are the three real roots of the cubic 
polynomial 4i^ — g2t — g^. The functions of the second type are periodic, with minimal pe- 
riod 2K{m), where K is complete elliptic integral of the first kind. The velocity of the traveling 
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waves originated by these functions is Vm = ~(1 + Tn{m — 1)). The first family of [41] consists 
of the traveling waves of the second type. When m — )• we obtain 

Ht) = ^ + ^csc{s + cf, 

and, when m — )• 1, we find 

3 13 

k{t) = 1 - -coth(s + c)^ k{t) = "2 + 2^<^ch{s + cf, 

which coincide with the solutions (67) and (69) of [38]. Next we consider the third-order equation 
k'" + kk' = 0. (6) 
Again, integrating twice, we find 

k'^ = -^k"^ - 12g2k - lUgs. (7) 
Each function satisfying (7) generates the traveling wave solution 
u{s,t) = k{s - 13252*) 

of the fifth-order Kaup-Kupershmidt equation. As in the previous case, the solutions of (7) can 
be expressed in terms of Weiertsrass p- functions and Jacobi elliptic functions: if A (5^2, 53) = 
-gl + 27gi > 0, then 

k{s) = — 12p(s + c), s £ {2nuji — c, (2n + l)iOi — c), n G Z, 

and, if A{g2,g3) < 0, we obtain 

k{s) = 4(1 + m) — 12ns(,s + c|m)^, k{s) = 4(1 — 2m) + 12mcn(s + c|m)^. 

The velocity of the traveling waves originated by these functions is Vm = — 176(1 + m{m — 1)). 
The second family of [41] consists of the traveling waves of the second type. When m — )• we 
obtain 

k{t) = 4 - 12csc(s + c)^ 
and, when m — )• 1, we find 

k{t) = 8 - 12coth(s + c)^ k{t) = -4 + 12sech(s + c)^ 
which coincide with the solutions (68) and (70) of [38]. 



3 Motion of curves in projective plane 

3.1 Curves in projective plane and their adapted frames 

Consider a smooth parameterized curve 7 : / RF^, defined on some open interval / C M 
and let G : I — )• M'^ \ {(0,0,0)} be any lift of 7. We say that 7(t) is an inflection point if 
Span{G{t),G'{t),G"{t)) has dimension < 2. From now on we will consider only curves without 
points of inflection. Then, 



r = Det{G,G',G")-^/^G 
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is the unique lift such that 

Det(r,r',r") = 1. (8) 

Differentiating (8) we see that there exist smooth functions a, 6 : I — t- M such that 
r'" = aT + bV. 

The projective speed v and the projective arc-element a are defined by 
V = {a - b' /2)^^'\ a = vdt. 

The primitives s : / — t- M of the projective arc-element are the projective parameters and the 
zeroes of a are the sextatic points. A curve without inflection or sextatic points is said to be 
generic. Obviously, every generic curves can be parameterized by the projective parameter. 

Remark 2. For every point t £ I, there is a unique non-degenerate conic Ct, the osculating 
conic, having fourth-order analytic contact with 7 at j{t). The osculating conic is defined by 
the equation — 2xqX2 = with respect to the homogenous coordinates of the projective frame 

"PI "P'l "P"! ^(^)"Pl 
J- i |j, i \t —i^ \t 

and, identifying MP^ with the space of plane conies, we define the osculating curve by 
C : tel -^C\te MP^ 

The sextatic points are critical points of the osculating curve. The assumption on the non- 
existence of inflection and sextatic points is rather strong from a global viewpoint. For instance, 
any simple closed curve in MP^ possesses flex or sextatic points and a simple convex curve has 
at least six sextatic points [35, 36]. 

The canonical projective frame field [3] along a generic curve 7 is the SL{3, ]R)-valued map 
defined by 

v' , 1 \ v' 1 ,, 

The canonical frame is invariant with respect to changes of the parameter and projective trans- 
formations. Furthermore, it satisfies the projective Frenet system 




F' = F-\l Q -k\v. (9) 
The function A; : / — )• M is the projective curvature, whose explicit expression [30] is 

where s is a projective parameter function and S is the Schwarzian derivative 
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Figure 1. The sextatic points and the osculating conic C(7r/4). 



Figure 2. The projective frame and the osculating conic. 



A 














u 











VJ 



Li 



Figure 3. The speed and the projective curvature. 

Remark 3. The construction of the canonical frame involves only algebraic manipulations and 
differentiations. So, it can be implemented in any software of symbolic calculus. In addition, 
the canonical frame can be constructed using the "invariantization" method of Fels-Olver [9]. 
In other words, there is a S'L(3, M)-equivariant map d : J^(M,MP2) 5L(3,M), defined on the 
fifth-order jet space of generic curves such that ^ o is the projective frame along 7, for 

every non-degenerate 7. 



Example 1. The convex simple curve 

- cos(t) , 



7 



t£R^ [(cos(t),sin(t),e~i~)] G 



has exactly six sextatic points, attained at 
n =0, 1.0412803807424216, 



vr, 



rs ^ 4.173209292276453, 



T3 ^ 2.109976014903134, 
tq ^ 5.241904926437165. 



Figs. 1 and 2 reproduce the curve, the sextatic points, the osculating conic and the projective 
frame at t = 7r/4. Fig. 3 reproduces the projective speed and the projective curvature. The 
speed vanishes at the sextatic points and the curvature becomes infinite at these points. 
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Figure 4. The projective curvature. 



Figure 5. The spherical lift of the corresponding curve. 



Remark 4. The curve is uniquely determined, up to projective congruences, by the speed 
and the curvature. If we assign smooth functions v > and k, the Frenet system (9) can be 
integrated with standard numerical routines (see Appendix A. 2). For instance, taking v = 1 
and choosing the "anomalous" 1-soliton solution of the /Ci-equation [17, 31] 



2m? (l + 2 cosh(m(s — niH)^ 



m 



0.8, 



2{2 + cosh{m{s-mH)Y) 

as projective curvature (see Fig. 4), the numerical solution of the linear system (9) gives rise to 
the curve whose spherical lift is reproduced in Fig. 5. 



3.2 The equation of a motion of curves 

A motion is a smooth one-parameter family 7(5, t) of projective curves such that 

7[i] : s G M ^ 7(s,t) G MP^ 
is generic and parameterized by the projective parameter, for every t £ I. Denoting by F^^-^ : 



SL{3, 



and k 



[t] 



the projective frame and the projective curvature of 7[t] we 



consider the projective frame and the projective curvature of the motion, defined by 
T: (s,t) G M X / ^ F[t](s) G S'L(3,M) 

and 

k: (s,t) G M X I ^kit]{s) G M. 
The projective frame satisfies 

T'^dT = /C(s, t)ds + $(s, t)dt, 
where 







— K 













— K 






Co 


1 








(10) 



(11) 



The coefficient 0q is said to be the normal velocity of the motion and it will be denoted by v. 
Proposition 2. The curvature of a motion of projective curves with normal velocity v satisfies 

dtK = Q{v,K) + XKs, (12) 

where X is a real constant, the internal parameter. Conversely, if k is a solution of (12) then 
there is a motion 7 with normal speed v, internal parameter A and curvature k. Moreover, 7 is 
unique up to projective transformations. 
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Proof. Differentiating (10) we obtain 

ds<^ - dtIC + [/C, = 0, (13) 
which imphes 

{4)s - KCPI + </.2 = 0, 

(,/.O),-2</,[]-0l + ^(0?-,/.i) = O, 

{(PDs + 0[] - -<Pl=0, 

(</)!), + + - cPl) -cPl = 0, 

{cf^Ds -4 = 0, 

(0?), + + k4 + 20} = 0, 

l)s - 3kcI)\ + 02 + ^1 _ ^ 0, (14) 



and 



dtK + (0i), - 0j + k(0O + 20}) + 0^ = 0. (15) 
If we set f = 0Q, then (14) gives 



18 1 \ fl 8 o\ /I 5 \ 5 1 

1^ + 7;KKs + TT^Ss + + -K 1^5+ ;; + T^Kg V2s + -KV3s + Tt:V5s, 



3 9 9 / V2 9y V6 6 / 9 18 



a X 1 . 2 ^^ /^l 4 2\ A 7 \ 4 1 

o = ^~qJ [kv3s + 4:K Vs)ds - \-K2s + -k ) ^ ~ ( 2 Is'^V '^^ ~ J)'^^^'' ~ is'"'^'" 

6? = -Xk + ^ y ('^'^^3s + 4K^fs)ds + + ^'^^ + ^'^s + ^'^s + ^/^^as + ^'^4s^ V (16) 



and 



'5 55 11 \ 4, 2 ^ /I 25 \ 11 1 

+ ( + Y^'^'^s + y'*3s I -(«^ + K2s)V2s+ I - + — I -^38+ yKt;4s+ ^^^^65, 



a 2 1 



.2 X 1 . 2 2 1 A A 7 \ 4 1 

-*1 = A / KU-js + 4k Vs)dS — -K H K2s f + Us Vs HV2s '^45, 

1 9 Jo V9 9 ; V2 18 ; 9 18 ' 



X 1 r / . 2 N , 2 2 \ 7 8 1 

^2 = A - - y (k-USs + 4k + l^-K + -K23 jV+ -l^sVs + g'«f^2s + -1^4^ 



(17) 



11 1 



a X r / . 2 ^ , Z^-, 4 3 1 82 1 \ 

02 = -Ak+-/ (kUSs + 4k f^jdS + ( 1 + -K - -K^ + + KK2s + gK4A 

/ 5 55 11 \ 4,2 N A 25 \ 

+ ("6''+ IS'"'"^ TS"^" jVs+^[n +H2s)v2s - ( g - Y^'^s ) 

where A is a real constant. From (16) and (17) we deduce that (15) is satisfied if and only if 
5jK = Q{v, k) + Xks- 

Conversely, if k is a solution of (12) and if we define 0*-, tC and <I> as in (11), (16) and (17), 
then K and ^> satisfy (13). Using Frobenius theorem we deduce the existence of a smooth map 

F : Rx I ^ SL{3,R) 

such that F-^dT = fCds + ^dt. The map is unique up to left multiplication by an element 
of 5L(3,M). Setting j{s,t) = [Fo(s,t)] we have a motion of projective curves with curvature k, 
normal velocity v, internal parameter A and projective frame J-. This yields the required result. 
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3.3 Local motions 

From the proof of Proposition 2 we see that the $-matrix of a motion of projective curves can 
be written as 

^ = ^{v,k) + XIC{k), 
where the coefficients of ^{v, k) are the integro-differential operators 




Bearing in mind the definition (1) of the hnear subspace -P[u], we deduce the existence of Hnear 
operators A^* : P[u] — )• J[u] such that 

for every p G -P[u]. This imphes the following corollary. 

Corollary 1. //p belongs to P[u\ and if k, is solution of the evolution equation 
dtK = 5(p)|j4^) + Xks 

then, there is a motion 7, uniquely defined up to projective transformations, with curvature k 
and normal velocity p|j^(K)- Motions of this type are said to be local. 

Remark 5. Local motions are the integral curves of local vector fields on the infinite-dimensional 
space V of unit-speed generic curves of MP'^. More precisely, if we take any p G -P[u] and any real 
constant A then there is a unique vector field ^p,A on V whose integral curve through 7^] S V 
is the local motion 7 such that: 
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• 7(s>0) = 7[o](s); 

• its curvature k is the solution of the Cauchy problem 

+ 5(p)|j^(«) + Xks = 0, k(s,0) = A;[o](s); 

• its normal speed is p|j^(K)- 

Definition 3. We say that -'^^p^A is the local vector field with potential p and spectral parame- 
ter A. The dynamics of a local vector field is governed by the induced evolution equation 

t^t + S{p)\j^(^^) + Xks = 0. 

From these observations and using Proposition 1 we have the following result. 

Theorem 1. For every n G N the local vector field X„^^Xn defined by the polynomial differential 
function On G P[u] and by the spectral parameter An induces the n-th equation of the Kaup- 
Kupershmidt hierarchy. 



4 Congruence motions 

Consider a local dynamics with potential p and internal parameter A. A curve 7 which evolves 
without changing its shape (by projective transformations) is said to be a congruence curve of 
the flow. Denote by k the curvature of 7 and by k{s, t) the curvature of the evolution 7(5, t) 
of 7(5). If k is non constant, then 

k{s, t) = k{s + vt), 

for some constant v. So, k is a traveling wave solution of the induced evolution equation and k 
satisfies the ordinary differential equation 

e(p|,4„),n) + (A + i;)n, = 0. (18) 

Unit-speed generic curves whose curvature satisfies (4) or (6) are examples of congruence curves 
of the first flow of the hierarchy. On the other hand, (4) is the Euler-Lagrange equation of the 
invariant functional defined by the integral of the projective arc-element a [3, 28]. This implies 
the following corollary. 

Corollary 2. Every critical curve of the functional 




is a congruence curve of the first flow of the Kaup-Kupershmidt hierarchy. 
The projective frame F{s, t) satisfies 

F'^dF = K{s + vt)dt + ^{s + vt)ds, (19) 

where the s[(3, M)-valued functions K and ^> are defined as in (11), (16) and (17), with normal 
speed p| • We define the Hamiltonian by 



H = ^-vK : M^5[(3,M). 
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The integrability condition of (19) is the Lax equation 

H' = [H,K] 
which implies the conservation law 

F ■ H ■ F-^ = 

where F is the projective frame of 7 and ^ is a fixed element of s[(3,M), the momentum of the 
congruence curve 7. In particular, H and ^ have the same spectrum. From now on we assume 
that F(0) = Id3x3. 

Proposition 3. The motion of a congruence curve 7 is given by 

jis,t) = ExpitO-7{s + vt). (20) 
Proof. Define 7(5, t) as in (20) and set 

F{s, t) = Exp(t^) • F{s + vt) V (s, t) G M X /. 

Since F is a lift of j{s,t), it suffices to prove that F satisfies (19). From the definition we 
deduce 

F-^dsF\(,^t)=K{s + vt) 

and 

F9tF|(,,t) = {F{s + vt)-^ ■ Exp(-tO) • (Exp(tO^ • F{s + vt) + vExp{tO ■ S^FU+^t) 
= F{s + vty^ F{s + vt) + vK{s + vt) = H{s + vt) + vK{s + vt) 
= ^{s + vt). 

This implies the required result. ■ 

We now prove the following proposition. 

Proposition 4. If k is a non-constant real-analytic solution of (18) and if ^ has three distinct 
eigenvalues then the corresponding congruence curve can he found by quadratures. 

Proof. It suffices to show that the solution of the linear system 

Fs = F-k, F(0) = Idsxs 

can be constructed from k by algebraic manipulations, differentiations and integrations of func- 
tions involving k and its derivatives. This can be shown with the following reasoning: 

Step I. The Hamiltonian H can be directly constructed from the prolongation j{k) of the 
curvature and from the potential p. Then, we compute the momentum ^ = H{0) and its 
eigenvalues tq, ti and T2. For each tj we choose row vectors H^^ and H^^ of H such that 

Sj = {'H^'-r,ej,)xCw--Tje,,)^0, 

where (eo, ei, £2) is the standard basis of C'^. 
Step II. Next, we define the S-matrix 

5= (^1,^2,53) : /^0[(3,C). 



Motions of Curves in the Projective Plane Inducing the Kaup-Kupershmidt Hierarchy 



15 



This is a real-analytic map which can be computed in terms of k and its derivatives. Subse- 
quently, we define the S-matrix 

E = F-S: M^0[(3,C). (21) 

The columns are eigenvectors of the momentum ^, with eigenvalue Tj, for each j = 0, 1, 2 

and every s. In particular, T,j has constant direction. Hence, there exist real-analytic complex 
valued functions rj such that 

S;. = r,S„ j = 0,l,2. (22) 

Differentiating S = F • 5 and using the structure equations satisfied by F we deduce 

S' + K-S = S-A{ro,ri,r2), 

where A(ro,ri,r2) is the diagonal matrix with elements rg, ri and r2. This shows that rp, ri 
and r2 can be computed in terms of k and its derivatives. 
Step III. We compute the integrating factors 

Pj{s) = Exp (^j^ rj{u)du^ , j = 0, 1, 2. 

Equation (22) implies 

^j = PjCj, i = 0,1,2, 
where Cq, Ci and C3 are constant vectors of C'^. We then have 

S = C- A(po,Pi,P2), (23) 
where C is a fixed element of GL(3,C). Substituting (21) into (23) we obtain 

F = M(0)-^ • M, 
where the M-matrix is defined by 

M-S = A{po,pi,P2). 

All the steps involve only linear algebra manipulations, differentiations and the quadratures of 
the functions tq, ri and r2, as claimed. ■ 

Example 2. We illustrate the integration of the congruence curves with projective curvature 

km{s) ~ ~ ~'~ + c|m)^) , m G (0, 1). 

Computing the iZ-matrix we obtain 

/h\ hi hf 
iJ= -2h\ hi 
V9 hi 

where the coefficients are given by 
3 

h\ = - (3dn(s|m)^ + m — 2) , h\ = Si — 9mcn(s|m)dn(s|m)sn(s|m), 
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9 

/if = - (-3dn(s|m)'' - 2(m - 2)dn(s|m)^ + m^) , 
,2 



hi = 9mcn(s|m)dn(s|m)sn(s|m) + 9. 

Next we compute the momentum and we get 

/|(m + l) 9 lim-lf' 

^= -3(m + l) 9 

V 9 |(m + l 

The discriminant of its characteristic polynomial is 

6 = —31 + m(6 + m(7 + (—6 + m)m)). 

Prom now on we assume 5 ^ Q. If 5 < the momentum has one real eigenvalue and two 
complex conjugate eigenvalues, otherwise the momentum has three distinct real eigenvalues. 
The eigenvalues are: 

^(m + n2Y'^ + 2^3 y/2{l - iy/2,){ni + na)^/^ + 2(1 + zV3)ra3 

^(ni + n2)V3 ' 2^(ni + n2)V3 

^(1 + iV3)(ni + n2)=^/^ + 2(1 - iV3)n3 
2^(m + n2)V3 



7-2 



where 

m = 3 (\/-3(-31 + m(6 + m(7 + (-6 + m)m))) - 9^ , 
n2 = (2 — m){m + l)(2m — 1), = {{m — l)m + 1). 

We set 

Sj = CH^ - Tjt2) X {'H^ - T,e3), j = 0, 1, 2, 
so that 

5j = - ^ (9dn(s|m)2 + 3(m - 2) - 2rj) (9dn(s|m)2 + 3(m - 2) + Tj) , 

S'j = 81 (mcn(s|m)dn(s|m)sn(s|m) + 1) , Sj=9 (9dn(s|m)^ + 3(m — 2)(/^ + r^) 

and 

9 — 9mcn(s|m)dn(s|m)sn(s|m) 
" ^ ^ ^ " 9dn(s|m)2 + 3(m - 2) + ' 

The quadratures can be carried out in terms of elliptic integrals of the third kind and we obtain 
pj = V2y^9dn(s|m)2 + 3(m - 2) + tj ■ e (3(-+i)+-i) ^ (24) 



where 



n(C;.^|m)= / " 
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is the incomplete integral of the third kind and am(u|m) is the amplitude of the Jacobi elliptic 
functions. Note that on the right hand side of (24) we have a smooth real branch of a multi- 
valued analytic function. The first column vector of is the transpose of 



, (to - n)(To - T2) ' (ri - To)(ri - T2) ' (t2 - To){t2 - Ti) 
and the coefficients M- of M{0)^^ are 



M, 



^ _ -9(m + 1)2 + 3(m + l)ri + 2Tf 
2y/6{m + 1) + 2Tj 



^2 ^ 81 ^3 _ 9V3(m + l)+r, 



3 



y/6{m + 1) + 2Tj ' ^ V2 



Then, the homogeneous components of a congruence curve with projective curvature km are 

and the evolution of the congruence curve is given by 
2 



{s, t) = Y, ^Mt^is^Hs - (1 + m(m - l))t), j = 0, 1, 2. 



fc=0 



Remark 6. These curves have a spiral behavior and a gnomonic growth (i.e. made of succes- 
sive self-congruent parts). Fig. 6 reproduces the spherical lifts of the congruence curves with 
parameter m = 0.3 and spectra 

ai = (-6.00233 + 6.06928i, -6.00233 - 6.06928i, 12.0047), 
CJ2 = (-9.94554, -9.94554, 19.8911), 

respectively. Fig. 7 reproduces the spherical lifts of the trajectories ■j(—,t) of the motion of 
a congruence curves with parameter m = 0.7, spectrum 

a = (-5.45852 + 6.40263i, 10.917, -5.45852 - 6.40263i) 

and t = 0, 0.5, 1, 1.5 respectively. 
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Figure 7. Trajectories 7(— ,t) of the congruence curve with parameter m = 0.3, spectrum a for t — 
0, 0.5, 1 and 1.5 respectively. 

A Appendix 

A.l Code to compute the equations of the Kaup— Kupershmidt hierarchy 

Jl[h_,y_]:=D[h,{s,3}]; 

J2[h_, v_]:=2 * D[v * Integrate[/i, s],{s, 2}]; 

J3[h_,v_]:=8 * 2 * Integrate + 3{v * D[h,s] + D[v * h,s]); 

J4[h_, v_]:=2 * Integrate[i; * D[h, {s, 2}] + 4 * v^2 *h,s]; 

J[h_, v_]:=Expand[FullSimplify[Jl[/i, + J2[/i, t;] + J3[/i, w] + J4[/i, u]]]; 

V[li.,Y.]:=D[h,{s,3}]+v*D[h,s]+D[v*h,s]; 

ii'[0][v_]:=l; H[l][v_]:=D[v, {s, 2}] + 4 * v^2; H[n_][y_]:=J[V[H[n - 2][v],v],v]; 

f) [n_] : =Expand [H [n] [u[s,t]]]; 

q[n_]:=Expand[Integrate[P[f)[n], t]], s]]; 

p [n_]:=Expand [Integrate [i?[n][e * u[s, t]] * u[s, t],{e, 0, 1]] 

KK[n_]:=Expand[i:'[u[s, t],t] + V[H[n][u[s,t]],u[s,t]]]; 

A. 2 Code to solve numerically the projective Prenet system 

Step I: define the speed, the curvature and domain of definition 
m:=0.8; t:=0; v[s.]:=l; a:= - 20; b:=20; 

,r 1 3m^(l+2Cosh[m*(g-m^4*t)]) . 

'*^P-J.— 2(2+Cosh[m*(s-m^4*i)])2 ' 

Step II: the routine to integrate the hnear system 

sol[l]:=NDSolve[{x'[t] == y[t],x[0] == l,y'[t] == -k[t] * x[t] + z[t],y[0] == 0, 
z'[t] == x[t] -k[t]* y[t],z[0] == 0}, {x, y, z}, {t, a, b}]; 

sol[2]:=NDSolve[{x'[t] == y[t],x[0] == 0,y'[t] == -k[t] * x[t] + z[t],y[0] == 1, 
z'[t] == x[t] - k[t] * y[t],z[0] == 0}, {x, y, z}, {t, a, b}]; 

sol[3]:=NDSolve[{x'[t] ==y[t],x[0] == 0,y'[t] == -k[t] * x[t] + z[t],y[0] ==0, 
z'[t] == x[t] - k[t] * y[t],z[0] == 1}, {x, y, z}, {t, a, b}]; 
5'[l][t_]:=Evaluate[{x[t],y[t],z[t]}/.sol[l]]; 
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5 [2] [t _] : =Evaluate [{x [t] , y [t] , z [t] } / . sol [2] ] ; 
S [3] [t _] : =Evaluate [{x [t] ,y[t], z[t]}/ . sol [3] ] ; 
r[t_]:={5[l][i][[l]][[l]],S[2][t][[l]][[l]],5[3][t][[l]][[l]]}; 
7[t-1:= , ^ Fftl; 
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